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ABSTRACT 

Context. XSS J 1227.0-4859 is a peculiar, hard X-ray source recently positionally associated to the Fermi- LAT source 1FGL J1227.9- 
4852/2FGL J1227. 7-4853. Multi-wavelength observations have added information on this source, indicating a low-luminosity low- 
mass X-ray binary (LMXB), but its nature is still unclear. 

Aims. To progress in our understanding, we present new X-ray data from a monitoring campaign performed in 2011 with the 
XMM-Newton, RXTE, and Swift satellites and combine them with new gamma-ray data from the Fermi and AGILE satellites. We 
complement the study with simultaneous near-UV photometry from XMM-Newton and with previous UV/optical and near-IR data. 
Methods. We analysed the temporal characteristics in the X-rays, near-UV, and gamma rays and studied the broad-band spectral 
energy distribution from radio to gamma rays. 

Results. The X-ray history of XSS J1227 over 7 yr shows a persistent and rather stable low-luminosity (6 x 10 33 d 2 kpc erg s _1 ) source, 
with flares and dips being peculiar and permanent characteristics. The associated Fermi- LAT source 2FGLJ1227. 7-4853 is also stable 
over an overlapping period of 4.7 yr. Searches for X-ray fast pulsations down to msec give upper limits to pulse fractional amplitudes 
of 15 - 25% that do not rule out a fast spinning pulsar. The combined UV/optical/near-IR spectrum reveals a hot component at ~ 13 kK 
and a cool one at ~4.6 kK. The latter would suggest a late-type K2-K5 companion star, a distance range of 1.4-3. 6 kpc, and an orbital 
period of 7-9 h. A near-UV variability (>6 h) also suggests a longer orbital period than previously estimated. 

Conclusions. The analysis shows that the X-ray and UV/optical/near-IR emissions are more compatible with an accretion-powered 
compact object than with a rotational powered pulsar. The X-ray to UV bolometric luminosity ratio could be consistent with a binary 
hosting a neutron star, but the uncertainties in the radio data may also allow an LMXB black hole with a compact jet. In this case, it 
would be the first associated with a high-energy gamma-ray source. 

Key words, accretion, accretion disks - X-rays: binaries - gamma rays: stars - binaries: close 


* Based on observations obtained with XMM-Newton , an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA, with Swift , a NASA science mis- 
sion with Italian participation, with Rossi-XTE, a NASA science mis- 
sion, and with Fermi a NASA mission with contributions from France, 
Germany, Italy, Japan, Sweden, and USA, and with AGILE, an Italian 
Space Agency mission with participation of the Italian Institute of 
Astrophysics and the Italian Institute of Nuclear Physics. 


1. Introduction 

XSS J12270-4859 (henceforth XSSJ1227) is a peculiar and 
enigmatic hard X-ray source. An early proposal that it could 
be a cataclysmic variable (CV) hosting a magnetic white 
dwarf (Masetti et al. 2006; Butters et al. 2008) was dis- 
regarded with independent observations by Pretorius (2009), 
Saitou et al. (2009), and de Martino et al. (2010, hereafter 
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dMIO). These observations instead suggested a low-mass X-ray 
binary (LMXB) nature. A putative 4.32 h orbital period was 
also claimed from optical and possibly near 1R photometry 
by dMIO, but not confirmed by further near-IR (nIR) data 
(Saitou et al. 2011). In dMIO we discovered this source to 
be positionally associated to the bright unidentified gamma-ray 
Fermi- LAT source, 1FGLJ1227. 9-4852, detected up to lOGeV. 
1FGLJ1227. 9-4852 is also reported in the 2nd Fenni-LAI cat- 
alogue as 2FGLJ1227. 7-4853 (hereafter 2FGLJ1227; Nolan 
et al. 2012). Based on this association, we found that the 
gamma-ray emission is a non-negligible fraction of the X-ray 
luminosity (L 0 .i-iooGev/Eo.2-tookeV ~ 0.8) and if the X-ray and 
gamma-ray emissions are linked, the peak energy should be be- 
tween l-100MeV. 

Due to its peculiar flaring and dipping behaviour, XSS J1227 
was proposed to share some similarities, though with rather dif- 
ferent timescales and energetics, with type-II burst sources and 
to represent an unusual low-luminosity (Lx ~ 6 x 10 33 erg s 1 at 
1 kpc) LMXB (dMIO). 

The nature of the compact object has been further inves- 
tigated in other wavelength domains (Saitou et al. 2011; Hill 
et al. 2011). From radio searches conducted by Hill et al. (2011), 
three sources within the Fermi- LAT 99% error box were de- 
tected, with only one having an obvious X-ray counterpart. 
This faint radio source is located at the position of XSS J1227 
with flux densities of 0.18 mJy at 5.5 GHz and of 0.14 mJy 
at 9 GHz. It was not detected at lower frequencies (640 MHz, 
240 MHz, and at 1400 MHz) in follow-up GMRT observations. 
For a power-law spectrum S v ~ v“, the radio data gave a fairly 
unconstrained power-law index a = -0.5 + 0.6. Also, Hill 
et al. (2011) did not detect fast radio pulses down to millisec. 
(msec) in subsequent Parkes observations. Nonetheless, the as- 
sociation with the Fermi- LAT source favoured XSS J 1227 as 
a msec pulsar (MSP) binary (Hill et al. 2011), similar to the 
first discovered rotational-powered MSP in a quiescent LMXB 
PSR J1023+0038 (Tam et al. 2010). The number of MSP bina- 
ries detected in gamma-rays by Fenni-LAI has rapidly increased 
in the last year with the newly discovered systems PSR J 1231- 
1411, PSR J0614-3329, and PSRI2214+3000 (Ransom et al. 
2011), PSR 11816+4510 (Kaplan et al. 20 1 2), PSR JO 1 0 1 -6422, 
PSR J1514-4946, and PSR J 1902-5 105 (Kerr et al. 2012). A 
few are found to host so-called black widow pulsars such as 
PSR J205 1 -0827(Wu et al. 2012) and J2339+0533 (Romani & 
Shaw 2011; Kong et al. 2012). A number of them are detected 
as faint X-ray sources with X-ray luminosities typical of radio 
MSPs (e.g. Ransom et al. 2011; Wu et al. 2012). Hence, it is 
expected that the Fermi satellite will discover most of the local 
population of MSPs, thus allowing constraint of the emission 
mechanisms and the binary evolution of their progenitors. 

On the other hand, from the detection of simultaneous X-ray 
and nIR flares and the similarity of broad-band spectral energy 
distribution (SED) at low energies with black hole (BH) can- 
didates such as GRS 1915+105, XTEJ1118+480 and CygX-3, 
Saitou et al. (2011) proposed that XSSJ1227 is reminiscent 
of microquasars with a synchrotron jet. However, broad-band 
SEDs, if due to synchrotron cooling in optically thin regime, 
are expected to show a break at ~500keV, making them un- 
detectable in the high-energy gamma-rays. The recent Fermi 
and AGILE detections of high-mass X-ray binaries (HMXBs; 
Abdo et al. 2009c; Tavani et al. 2009b; Ackermann et al. 2012; 
Hadasch et al. 2012), with some of them still thought to host 
neutron stars (NS) or BH (Hadasch et al. 2012; Papitto et al. 
2012; Torres et al. 2012) may open the possibility of detecting 


faint gamma-ray sources in different binary types (Hadasch et al. 

2012). 

To further advance our understanding of this intriguing ob- 
ject, we started a X-ray monitoring campaign in 2011 to infer 
whether this source undergoes changes in luminosity and to de- 
tect possible fast (msec) pulsations in the X-ray band. We were 
also granted AGILE/GRID data acquired since the beginning 
of operations in 2007 to detect for the first time this source in 
the soft gamma rays down to 50 MeV to obtain a broad-band 
gamma-ray coverage combining the GRID and LAT data. With 
these aims, we present the X-ray observations acquired along 
one year with Swift and Rossi-XTE (RXTE) and the fast timing 
X-ray data acquired with XMM-Newton. In Sect. 2 we report the 
observations and data reduction; in Sect. 3 we provide analysis 
of the X-ray light curves, including the long-term comparison 
with the gamma-ray Fermi- LAT curve and the search for fast 
msec pulses; in Sect. 4 we examine the spectral properties over a 
wide energy range, and in Sect. 5, we discuss the possible nature 
of this object. 

2. Observations and data reduction 

The new observations acquired with XMM-Newton, RXTE, and 
Swift as well as the time coverage of AGILE/GRID are sum- 
marised in Table 1 . 

2. 1. The XMM-Newton observation 

A 30 ks pointing with XMM-Newton (OBSID: 656780901) was 
carried out on Jan.l, 2011. In contrast to our past observation 
in 2009 (dMIO), the EPIC-pn camera (Striider et al. 2001) was 
operated in the timing read-out mode that reaches a resolution 
of 0.03 ms. The EPIC-MOS cameras (Turner et al. 2001) were 
instead operated in imaging full-window mode using the thin fil- 
ters. The optical monitor (OM; Mason et al. 2001) was operated 
in fast-window mode using the U (3000-3800 A) filter through- 
out the observation (see Table 1). The data were processed us- 
ing the standard reduction pipelines and analysed with the SAS 
10.0 package using the latest calibration files. The photon ar- 
rival times from EPIC cameras and OM were reported to the so- 
lar system barycentre using the nominal position of XSS J 1227 
(Masetti et al. 2006). In the timing-mode read-out used for the 
EPIC-pn exposure, the spatial information is lost because imag- 
ing is made only in one dimension and the data from a predefined 
area on one CCD chip are collapsed into a one-dimensional row 
to be read out at high speed. The EPIC-pn events from the source 
and background were extracted from RAWX = 24-52px and 
from RAWX = 3-15px, respectively. For the EPIC-MOS cam- 
eras, we instead extracted events using a circular region of 37" 
aperture radius centred on the source and using a background 
region located on the same CCD chip. In order to improve the 
signal-to-noise ratio (S/N), we filtered the data by selecting pat- 
tern pixel events up to double with zero-quality flag for the 
EPIC-pn data, and up to quadruple pixel events for the EPIC- 
MOS data. The average background level of the EPIC cam- 
eras was low during the whole observation. The background- 
subtracted OM-U light curve was obtained with a binning time 
of 10 s. 

2.2. The RXTE observations 

In order to infer the mid- to long-term behaviour of XSS J1227, 
we initiated a monitoring with RXTE (Bradt et al. 1993, 


A89, page 2 of 16 


D. de Martino et al.: Follow-ups of XSS J12270-4859 


Table 1. Summary of the observations of XSS J1227. 


Telescope 
Prog. ID 

Instrument 

Date 

UT (start) 

Exposure time (s) 

Net count rate (counts s *) 

XMM-Newton 4 

EPIC-pn 

2011-01-01 

07:06 

30045 

3.71 ±0.09 

65678 

EPIC-MOS 


06:47 

31022 

1.09 ±0.01 


OM-U 


06:56 

2699 

5.19 ± 0.04 




07:46 

2700 





09:07 

2700 





09:57 

2699 





10:48 

2700 





11:38 

2499 





12:25 

2499 





13:12 

2499 





13:59 

2498 





14:46 

2500 


RXTE 5 

PCA 

2011-01-03 

00:27 

2013 

1.72 ±0.06 

96309 


2011-01-10 

14:11 

2000 

1.83 ±0.06 



2011-01-16 

17:38 

1917 

1.94 ±0.05 



2011-01-23 

12:37 

2082 

2.21 ±0.07 



2011-01-31 

13:26 

1919 

2.03 ±0.07 



2011-02-07 

13:38 

2052 

0.96 ±0.06 



2011-02-14 

15:05 

1947 

1.84 ±0.06 



2011-02-21 

17:49 

1836 

0.97 ± 0.06 



2011-02-28 

12:51 

2476 

1.59 ±0.06 



2011-03-07 

10:58 

1977 

2.44 ± 0.07 



2011-03-14 

05:59 

2043 

2.05 ± 0.07 



2011-03-21 

12:32 

1962 

0.89 ±0.06 



2011-03-28 

21:22 

1874 

0.56 ±0.06 



2011-04-03 

18:31 

1794 

1.07 ±0.06 



2011-04-13 

09:03 

1919 

1.73 ±0.06 



2011-04-18 

14:18 

1770 

2.25 ± 0.07 



2011-04-25 

06:12 

1962 

1.90 ±0.06 



2011-05-03 

19:58 

1586 

1.05 ±0.08 



2011-05-10 

21:05 

1922 

1.82 ±0.07 



2011-05-16 

13:17 

1909 

1.03 ±0.06 



2011-05-23 

17:37 

2010 

1.63 ±0.06 



2011-05-30 

15:49 

2034 

0.91 ±0.06 



2011-06-06 

06:14 

1742 

2.06 ± 0.07 



2011-06-13 

02:34 

1807 

1.09 ±0.06 



2011-06-19 

04:14 

2043 

2.27 ± 0.07 



2011-06-27 

17:23 

2099 

1.83 ±0.07 



2011-07-04 

14:18 

1824 

1.34 ±0.06 



2011-07-11 

10:34 

1789 

1.06 ±0.07 



2011-07-18 

13:33 

1965 

1.51 ±0.07 



2011-07-25 

11:39 

1716 

1.51 ±0.06 



2011-08-01 

03:25 

1947 

2.06 ± 0.07 



2011-08-08 

06:32 

2102 

1.94 ±0.06 



2011-08-15 

06:00 

2004 

0.84 ±0.06 



2011-08-22 

08:38 

1756 

0.90 ±0.07 



2011-08-29 

17:41 

2132 

1.61 ±0.06 



2011-09-19 

21:22 

1859 

1.10 ±0.06 



2011-09-27 

04:51 

1656 

0.83 ±0.06 



2011-10-08 

12:01 

1731 

1.03 ±0.07 



2011-10-21 

20:53 

2025 

1.03 ±0.06 



2011-11-10 

18:09 

1933 

2.24 ±0.07 



2011-11-20 

14:50 

1860 

0.80 ±0.06 



2011-11-30 

18:27 

2094 

0.84 ±0.07 



2011-12-10 

00:03 

1787 

0.84 ±0.06 



2011-12-23 

02:12 

1947 

1.04 ±0.06 


Notes. (4) Total EPIC count rates in the 0.2-10keV range (5) Total RXTE/PCA count rates in the 2-9 keV range. (6) Total Swift/XKT count rates in 
the 0.3-10keV range. 171 AGILE/GRID Pointing and spinning mode data: the start and end time and total exposure are reported. 
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Table 1. continued. 


Telescope 
Prog. ID 

Instrument 

Date 

UT (start) 

Exposure time (s) 

Net count rate (counts s ') 

Swift 6 

XRT 

2011-03-23 

04:35 

672 

0.17 ±0.09 

41135 


2011-03-23 

06:11 

732 

0.26 ±0.11 



2011-03-23 

07:47 

701 

0.26 ±0.09 



2011-03-23 

09:24 

729 

0.28 ±0.12 



2011-03-23 

11:00 

729 

0.35 ±0.12 



2011-03-23 

12:36 

852 

0.27 ±0.09 



2011-03-23 

14:13 

792 

0.30 ±0.09 



2011-05-10 

16:32 

634 

0.67 ±0.15 



2011-09-15 

07:55 

854 

0.32 ±0.13 



2011-09-18 

09:48 

809 

0.07 ±0.06 



2011-09-19 

04:59 

1483 

0.39 ±0.15 



2011-09-19 

06:34 

859 

0.17 ±0.09 



2011-09-19 

08:12 

802 

0.22 ±0.10 



2011-09-19 

22:39 

392 

0.34 ±0.14 



2011-09-21 

02:02 

654 

0.31 ±0.13 



2011-09-21 

10:11 

242 

0.15 ±0.08 



2011-09-22 

21:16 

557 

0.28 ±0.12 



2011-09-25 

02:11 

1293 

0.39 ±0.15 



2011-09-26 

02:20 

557 

0.24 ±0.09 

AGILE 7 

GRID 

2007-10-01 

2009-10-31 

5.3 x 10 6 


70 


2009-11-15 

2011-05-15 

3.8 x 10 6 



Prog. ID: 96309) that started in Jan. 2011 and ended in 
Dec. 2011. A total of 44 RXTE/PCA pointings were performed 
with a cadence of about one week and typical exposure times 
of 2 ks (see Table 1). The RXTE/PCA standard data products for 
each observations were obtained from the RXTE Guest Observer 
facility. Barycentric corrections were applied to the background- 
subtracted light curves. XSS J1227 was typically found at a net 
count rate of ~0.3-0.5ctss -1 /PCU in the 9-20 keV. Due to the 
poor statistics above lOkeV, the timing analysis is restricted 
to the 2-9 keV range, also using the soft 2-4 keV and hard 
4-9 keV bands. The source spectra were instead analysed up to 
30 keV in each observation. 

2.3. The Swift observations 

The monitoring of XSS J1227 with Swift [X RT (Prog.ID:41 135) 
started in Mar. 2011 and ended in Sept. 20 1 1 . A total of 19 snap- 
shots were performed with typical exposure times of 800 s. A 
denser coverage was performed on Mar. 23 and on Sept. 19 (see 
Table 1). We used the Swift /X. RT data products generator at the 
University of Leicester (Evans et al. 2009) to build background- 
subtracted light curves in the 0.3- lOkeV, 0.3-1.5keV, and 
1.5-10keV bands and spectra in the 0.3-10keV range. We used 
photon counting-mode data only. Barycentric corrections were 
applied to the extracted light curves. 

2.4. The AGILE data 

Astrorivelatore Gamma ad Immagini LEggero (AGILE; Tavani 
et al. 2009a) consists of a large field of view (~60°) gamma-ray 
imager, GRID, sensitive in the energy range 30MeV-50GeV, 


and a co-axial hard X-ray detector (Super-AGILE) for imaging 
in the 18-60keV range (Feroci et al. 2007). Since beginning op- 
eration in Jul. 2007, AGILE has worked in pointing mode, but 
science operations were reconfigured following a malfunction 
of the rotation wheel in mid-Oct. 2009. Since then, the satel- 
lite has been operating in “spinning observing mode”, with the 
solar panels pointing at the Sun and the instrument axis sweep- 
ing the sky with an angular speed ~1 deg s -1 . After our associa- 
tion of XSS J1227 to the Fermi source (dM10), we were granted 
AGILE/GRID data in Cycle 3 (Prog.ID:70). We have enlarged 
the dataset including observations collected in pointing mode 
from Oct. 1, 2007 to Oct. 31, 2009, and the data collected with 
the satellite in spinning mode from Oct. 31, 2009 to May 15, 
2011; the total exposure times amount to 5.3 Ms and 3.8 Ms 
respectively. 

XSS J 1227 is too faint in the X-rays to be detected by 
the Super-AGILE instrument and therefore only the GRID 
dataset was analysed using the Build 20 release of the GRID 
analysis software available at the ASI-Science Data Centre 
(SDC). Well-reconstructed y-ray events were selected using 
the FM3.119 filter. Events collected during the passage in 
the South Atlantic Anomaly and Earth-albedo gamma rays 
coming from a circular region of radius 80° centred on the 
Earth were rejected. A maximum likelihood analysis was per- 
formed in the 100MeV-50GeV energy band, obtaining 2 cr up- 
per limits of 7 x 1CL 8 photons crcr 2 s~* with a total exposure 
coverage of 2 x 10 y cm 2 s for the pointing-mode period and 
1.2 x 1CL 7 photons cm~ 2 s _1 with a total exposure coverage 
of 1.5 x 10 9 cm 2 s for the spinning-mode period. For the to- 
tal period under analysis we obtained a 2cr upper limit of 
6 x 1CL 8 photons cnT 2 s _1 with a total exposure coverage of 
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~3.4 x 1 0 y cm 2 s. In order to investigate soft (E < 200 MeV) 
gamma-ray emission from XSSJ1227, we also integrated the 
GRID data in the non-standard 50-200 MeV energy band, ob- 
taining a 2cr upper limit of 14.8 x 10~ 8 photons cm -2 s~' with a 
total exposure coverage of 2 x 10 9 cm 2 s. 

2.5. The Fermi data 

The Fermi- L AT data were retrieved from the Fermi Science 
Support Centre. The dataset spans 44 months since the start 
of Fermi operations, i.e. from Aug. 4,2008 to Apr. 17,2012, 
amounting to a total exposure time of 52.6 Ms. We used Pass7 
photon data, reducing and analysing them using the Fermi 
Science Tools v9r27 package. We used the high-quality (dif- 
fuse) photon event class (EVENT CLASS = 2) and the Pass 7 
v6 Source (P7 V6 source) instrument response functions (IRFs). 
We excluded time periods when the region around 2FGL J1227 
was observed at a zenith angle greater than 105° to reduce con- 
tamination by Earth albedo gamma rays. Correction to the solar 
system baricentre was also applied. 

We produced a LAT aperture photometry light curve with 
a time bin of four days from a circular region centred on 
the source position and a radius of 1° to avoid contamination 
from close sources (2FGLJ121 8.8-4827, 2FGL J1231.3-51 12 
and 2FGLJ1207. 3-5055) (see Nolan et al. 2012). We performed 
this operation with the gtbin task and applied the exposure cor- 
rection using the gtexposure task. These tasks do not perform 
background subtraction, and hence no background correction 
was applied. 

The spectral fluxes in the five LAT bands (100-300 MeV, 
0.3-1 GeV, 1-3 GeV, 3-10GeV, and 10-100 GeV) were instead 
collected from the 2nd source catalogue that covers 24 months of 
operations and were obtained using maximum likelihood analy- 
sis (see details in Nolan et al. 2012). 2FGLJ1227 is reported 
at a flux of 3.34 + 0.23xl0~ n ergcm~ 2 s _1 in the 0.1-100 GeV 
range at a significance of 24.3 cr in accordance with the value of 
3.95 + 0.44 x 10 11 erg cm 2 s 1 reported in the 1st LAT source 
catalogue. The localization procedure used to construct the cat- 
alogue also provides spectral fits to all sources. In contrast to 
the 1 st catalogue, which used only a power-law function, differ- 
ent spectral shapes are used in the 2nd one. For 2FGLL1227, a 
power-law and a function, called “logParabola” are fitted to the 
spectrum. This representation, allows for a smoother decrease 
at high energy than a power-law exponential cutoff form. It re- 
duces to a simple power-law when the curvature parameter [5 = 0 
and is reported in the catalogue when the curvature is significant 
by above 4cr (Nolan et al. 2012). A simple power-law fit is re- 
ported with a spectral index of 2.3306. The spectral index of the 
logParabola is 2.085 ±0.087, f3 = 0.288 + 0.064 with a pivot en- 
ergy at 549.30 MeV (where the uncertainty on differential flux is 
minimal). The significance of the fit improvement with respect 
to the simple power-law is 5.6cr. 

3. Results 

3. 1 . X-ray persistent characteristics 

The background-subtracted X-ray light curves from EPIC-pn 
and MOS cameras were first extracted in 60s bins to evaluate 
the overall behaviour of XSSJ1227 during the XMM-Newton 
observation. The light curves from the EPIC-pn and MOS are 
similar, and we use the EPIC-pn data because of the higher S/N. 
In Fig. 1, the 0.3-10keV light curve with a binning time of 10 s 
is shown together with the temporal behaviour of the hardness 


ratio HR = [H-S/H+S], where S and H are the count rates in the 
0.3-2 keV and 2-10keV, respectively, with a binning time of 
60 s. XSS 11227 was found to be at a persistent (quiescent) level 
of ~4.35 ± 0.35 cts s~ 1 for most of the EPIC-pn exposure, except 
for the last 4 ks, where the count rate increases up to a factor 
of 6.4. While the general behaviour observed in 2011 is simi- 
lar to that of the Jan. 2009 XMM-Newton observation (dM10), 
the EPIC-pn timing mode data show that the persistent level is 
characterized by dips of variable length (see Sect. 3.3). 

The hardness ratios confirm the general trend observed in 
2009 (dM10), where a substantial hardening is detected dur- 
ing the postflare dips and there is no major spectral change 
during flares and quiescent dips. The intensity-versus-hardness 
ratio diagram (HID) in Fig. 2 is similar (see Fig. 4 of dM10). 
Here we define flares when the count rate >5.7 cts s - 1 (cyan in 
Fig. 2) and dips have a count rate <3.ctss~'. The higher S/N 
EPIC-pn data allow the dips to be better sampled. Therefore 
we further separate soft dips with HR < -0.6 (blue) and hard 
dips with HR > -0.2 (red) from “ordinary” neutral dips with 
-0.6 < HR< -0.2 (magenta). The hard and soft dips are defined 
as those below/above 2 cr the average HR value in quiescence. 
The HID confirms that only post-flare dips are hard, while dips 
during quiescence are mostly neutral with a marginal softening 
in those that are longer. 

Flares and dips in the X-ray light curve of XSS J1227 are 
found in all observations conducted in 201 1 by RXTE and Swift. 
In Fig. 3, we show the light curves of each RXTE pointing 
with a binning time of 16 s. The hardness ratios HRrxte = 
[H - M/H + M] in the 2-4 keV and the 4-9 keV bands were 
constructed with a binning time of 128 s. We define the qui- 
escent level as that observed on Dec. 23 (Table 1), and hence 
flares are defined with a count rate >1.8 cts s^ 1 and dips with a 
count rate <0.8 cts s~' . The average hardness ratio in quiescence 
is HRrxte = 0.38 ± 0.24, and we define hard/soft dips as those 
over/below 2cr this average value, respectively. Ordinary neutral 
dips are instead those found within this range. The RXTE HID 
(Fig. 4) shows that flares and ordinary and soft dips have been 
relatively well sampled, but not the hard dips. This is likely due 
to the higher energy coverage of RXTE/PCA with respect to 
XMM-Newton and Swift instruments. 

The Swift light curves (not shown) were extracted with 
a binning time of 25 s. Again, XSS LI 227 is found stable at 
~0.2-0.3 cts s _1 during most of Swift pointings (see Table 1). 
Despite the poor coverage, several flares are detected or partially 
detected, with one notably on Sept. 19 that reached peak inten- 
sity ~8 times the quiescent (0.26 cts s~') level. The flares are de- 
fined with a count rate >0.43 cts s -1 , while dips are defined with 
a count rate <0.2 cts s^ 1 . The hardness ratios were constructed 
with a binning time of 128 sin the 0.2-1. 5keV and 1.5-10keV 
bands. Soft/hard dips were defined as over/below 2cr the average 
value HR Sw(ff = 0.1 +0.1. The HID from Swift also reveals the 
tendency of hardening during post-flare dips (Fig. 5). 

We also inspected archival Swift unpublished observations 
on Sept. 15, 2005, Sept. 24, 2005 and Aug. 10, 2010 in similar 
fashion as the 2011 Swift pointings to construct a long-term his- 
tory of XSS J1227. 

The data acquired during the 2011 monitoring compared 
with the observations carried out in 2005 by Swift, in 2007 
by RXTE (dM10), in 2008 by Suzcikii (Saitou et al. 2009), 
in 2009 with XMM-Newton (dM10) and RXTE (Saitou et al. 
2011), and in 2010 by Swift allow to definitively assess that 
XSS J 1227 is a persistent and not a transient source. It is stable 
on a 7-yr timescale while undergoing repetitive flaring episodes. 
From the longer datasets we also estimate flare occurrence and 
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Fig. 1 . Bottom panel: net U band light curve with a 10 s binning time. Count rates are normalized to the average persistent quiescence value. A 
sinusoidal fit with a period of 6.4 h is also reported. Gaps are due to the OM fast acquisition mode. Middle panel: the EPIC-pn net light curve in 
the 0.3-10keV range with a 10 s binning time. The count rate is normalized to the average quiescent value. Flares are marked in cyan, dips with 
no spectral change in magenta, dips showing a softening are marked in blue, while hard dips are marked in red. Arrows mark ten long dips used 
to study the CCF and spectra. Top panel: the hardness ratio between the 0.3-2 keV and 2-10 keV bands with binning time of 60 s. (See electronic 
edition for colour version.) 


duration: Aq ule 2 24 ks and Apiares ~ 8.2 ks from Suzaku (2008); 
AQuie ~ 35 ks and Apiares ~ 3.6 ks from RXTE (2007) and 
AQ U ie = 27.1 ks and Apices S 3.1 ks from XMM-Newton (2009) 
and AQ U i e ~ 27.4 ks and Api^s ^ 3.4 ks from XMM-Newton 
(2011), giving lower limits of AQ U j e = 27-35 ks and Api ares ~ 
4ks, the latter grouped in multi-flares. These could be a signa- 
ture of a ^9-11 h occurrence. 


3.2. The near-UV light 

Figure 1 also reports the simultaneous light curve recorded with 
the OM in the U band filter in Jan. 2011, where the source also 
undergoes a flaring activity at the end of the XMM-Newton ob- 
servation. XSSJ1227 is at an average U band magnitude of 
16.6 mag during quiescence and reaches 15.1 mag at flares. 


The highest peak in the U light curve is ~3 times the persistent 
level, smaller than the X-ray variations. Despite the gaps in this 
light curve, it is apparent that dips occur also in the near-UV 
during quiescence, with the count rate dropping by a factor of 
~1.4. Hence, in contrast to what was claimed in dM10, it is 
now possible to identify the UV counterparts of the X-ray dips 
in XSS J1227. An enlargement of the light curves is shown in 
Fig. 6. 

Furthermore, a long-term variation in the UV was suspected 
to be present in the 2009 XMM-Newton OM data. A Scargle pe- 
riodogram of the quiescent U band time series shows a peak 
at 6.7 h (at 3.6cr c.l.). A formal sinusoidal fit gives a period of 
6.4 + 0.2 h (uncertainty at lcr c.l.) and full amplitude of 22 ± 1% 
(see also Fig. 1). The statistical significance of the inclusion of 
a sinusoid with respect to a constant function, evaluated with 
a F-test, is 10cr. Given that the length of the U band exposure 
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Fig. 2. HID between 0.3-2 keV and 2-10keV EPIC bands. Flares are 
marked in cyan, dips with no spectral change in magenta, soft dips in 
blue, and hard dips in red. Green points represent quiescence. (See elec- 
tronic version for colour figure.) 

during quiescence is 26 ks, the 6.4 h period should be regarded as 
a lower limit. We note here that, in contrast to what was claimed 
from the past XMM-Newton pointing (dMIO), no variations on 
timescales of hours are found in this EPIC-pn quiescent light 
curve. 

3.3. The dips in focus 

The X-ray dips represent a persistent characteristic of XSS J1227 
during quiescence. The time spent in the dips amounts to 7.78 ks, 
representing ~30% of quiescent period. A few (ten are marked 
in Fig. 1) longer dips have durations between 100 s and 700 s, 
but the one observed prior the flaring activity is the longest, with 
a duration of 1 100 s. Similarly, in the 2009 XMM-Newton point- 
ing, the longest dip is detected before the flaring activity. In both 
observations, these dips are soft. The source also undergoes sev- 
eral short dips of length 10-90 s, where the count rate does not 
drop to zero. These short dips are the ordinary neutral dips with 
no substantial spectral change (Fig. 2). 

The selected ten long dips are characterised by a less steep, 
though rapid, decay reaching almost zero counts in the longer 
ones and a faster recovery to the persistent level. The decay and 
rise times have lengths ranging between -40-90 s and appear to 
be correlated with the length of the dips, with slower decays/rises 
in the longer dips. The hardness ratios show a softening only in 
the longer dips. We also performed a cross correlation (CCF) 
between the hard (2-10keV) and soft (0.3-2 keV) light curves 
during the ten selected long dips. The CCF is shown in the right- 
hand lower panel of Fig. 7, where the soft band is taken as a ref- 
erence and its autocorrelation function (ACF) is also reported. 
The CCF does not reveal any lag between the two bands, imply- 
ing that the bulk of X-ray flux arises from the same region. 

Near-UV dips seen in quiescence occur almost simultane- 
ously with the X-ray dips, but the flux deficiency is ~40%, hence 
much lower than their X-ray counterparts (Fig. 6). The shape of 
UV dips is also different, being shallower with a smooth decay 
and rise. This is reflected in the count rate ratios between the two 
bands. Using the total X-ray range as reference, we computed 


its ACF and cross-correlated the U band light curve with the 
X-ray one using eight of the ten dips covered simultaneously 
(see Fig. 1). The CCF is shown in the right-hand top panel of 
Fig. 7 together with the ACF of the X-ray band. It is peaked at 
zero lags, but an excess at negative lags ~1 00-300 s is found. 
Both the UV deficiency and CCFs indicate that the UV-emitting 
region is more extended. 

3.4. The flaring states 

XSS J 1227 showed flares at the end of the XMM-Newton obser- 
vation where the X-ray count rate increased by a factor of ~6.4 
(Fig. 8). Similar to the 2009 XMM-Newton light curve, the flar- 
ing activity consists of multiple events. Here we identify at least 
four major events, lasting about 600 s, 640 s, 847 s, and >350 s 
(the latter covered only at the start). The higher S/N EPIC-pn 
timing mode data, while confirming the general exponential rise 
and rapid decay of major flares, also allow us to identify that 
they are highly structured and composed of multiple peaks with 
a typical timescale of 30-50 s. The major spectral changes oc- 
cur during the dips after these flares, as depicted in Fig. 1, where 
hardening is detected but not at flare rise and maximum. We 
cross-correlated the hard (2-10keV) with the soft (0.3-2keV, 
taken as reference) bands and computed the ACF in the soft band 
in the flare portions during the rise and maximum. The CCF and 
ACF are shown in Fig. 7, where no delay is inferred, confirming 
what was found in 2009. 

The flares, simultaneously detected in the U band, display 
peak intensities up to a factor of ~3. In the U band, only one 
post- flare dip can be recognized precisely after the third flare 
(Fig. 8). The CCF between U and X-ray (taken as reference) 
peaks at zero lags but is also asymmetric towards negative lags 
(-150-200 s). It is different from what was found in 2009, possi- 
bly because of the short coverage between the two bands in that 
observation. This could be due to the rather different structure of 
U flares with respect to the X-ray ones, where not all post-flare 
X-ray dips have their UV counterpart and the last UV flare rises 
faster than the X-ray one. 

3.5. Search for fast coherent signals 

The EPIC-pn timing mode data were also inspected to search for 
fast periodic or quasi-periodic signals in the range 0.5-10 keV. 
In order to assess the presence of a coherent signal in the EPIC- 
pn observation of XSS 11227, a Fourier analysis was first per- 
formed using the whole exposure, binning the time series at eight 
times the temporal resolution of the EPIC-pn timing mode data, 
so that fbin — 0.236 ms. Hence the search for periodic signals 
was performed in the frequency interval 3.3 x 1 0 5 — 2 1 1 4 Hz. 
We took into account the number of frequency bins searched 
to define a threshold power level to be exceeded in the pres- 
ence of a signal above the noise (see, e.g. van der Klis 1989). 
Because of the presence of red noise at low frequencies (Fig. 9), 
we used a smoothing window technique to evaluate the detec- 
tion level at each of the frequencies searched (Israel & Stella 
1996). No signal was significantly detected at a confidence level 
of 3cr. Given the maximum power detected in the power spec- 
trum, we evaluated an upper limit at 3cr confidence level on the 
fractional amplitude of pulsations of /U) L = 0.050. The effects of 
the observation length and binning, as well as the interaction be- 
tween a putative signal and noise, are accounted for as described 
by Vaughan et al. (1994). The power spectrum continuum be- 
low 100 Hz is shown in Fig. 9. Three flat-top noise components 


A89, page 7 of 16 



A&A 550, A89 (2013) 



Fig. 3. RXTE net light curves during the pointings in 201 1 with a binning time of 16 s. The count rate is normalized to the quiescent average value 
observed on Dec. 23. Flares and dips are shown in colours as in Fig. 1. (See electronic version for colour figure.) 


of width, Wi = (1.0*$ x 1(T 2 , W 2 = 0.19+°-°®, and W 3 = 
4.0*J q Hz, adequately fit the power spectrum (y 2 = 37/67). 

The upper limit on the pulsation amplitude is meaningful 
only if a putative pulsar is isolated or belongs to a wide bi- 
nary system (P or b » 7’ 0 bs). In dM10 we reported a possi- 
ble 4.32 + 0.01 h optical photometric periodicity, possibly asso- 
ciated with the binary orbital period. The XMM-Newton OM 
U band light curve indicates the presence of a longer period 
6.4 + 0.2 h (see also Sect. 4). If the binary period is as short 
as 4. 3—6.' 4 h, the secondary star mass would be in the range 
0.3-0. 7 M q (Smith & Dhillon 1998). We therefore need to con- 
sider shorter integration times in order to limit the leakage of 
power in close frequency bins due to the Doppler shift affecting 
the signal frequency. The optimal integration time to search for a 
signal coming from a binary system, without making any correc- 
tion of the photon arrival times for the unknown orbital motion, 
was evaluated by Johnston & Kulkarni (1991). Adopting as a 
lower limit the 4.3 h period and using the corresponding orbital 
parameters, one obtains 


( v Pi 

f Porb fl 
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( sin i \ 
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Fig. 4. HID between the HRrxte 2-4 keV and 4-9 keV bands as de- 
fined in the text. Colours are used as in Fig. 2. (See electronic version 
for colour figure.) 

( 1 ) 
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Fig. 5. HID between the 0.3-1.5 keV and 1.5- lOkeV SwiftfXRT bands. 
Colours are used as in previous figures. (See electronic version for 
colour figure.) 



Time(s) since HJD 2455556S.8 

Fig. 6. Enlargement of the normalized light curves in 201 1 showing the 
simultaneous occurrence of quiescent dips in the U band ( bottom panel) 
and X-rays ( middle panel) (0.3- 10 keV) with a binning time of 10 s, to- 
gether with the ratio of the U and X-ray band count rates (top). Vertical 
dashed (red) lines mark the dip start and end times. 


Assuming that a pulsar in this LMXB system has a frequency 
in the range displayed by the known accreting pulsars with such 
low mass companions (100-600 Hz), one obtains optimal inte- 
gration times spanning 150-395 s. We produced several power 
spectra over time intervals of length 124 and 248 s 1 2 . No sig- 
nificant detection within a 3 cr confidence level was obtained. 
However, this is not surprising however, because at the low count 

1 The lengths of the time intervals 124 s, 248 s, and 1983.7 s were cho- 
sen in order to give a number of bins equal to an integer power of 2 19 , 

2 20 , and 2 23 , respectively. This greatly decreases the CPU time needed 
to perform fast Fourier transforms (FFT). 


Fig. 7. CCFs (solid lines) during flares (left) and dips (right). The CCF 
between the hard and the soft X-ray bands, taking as reference the soft 
range (ACF), are shown in the bottom panels. The CCFs between the 
total X-ray and U bands, taking as reference the X-ray band (ACF), are 
shown in the top panels. The ACF of each reference band is plotted with 
a dotted line. 



Time(s) since HJD 24555562.8 

Fig. 8. Enlargement of the normalized light curves showing the flar- 
ing activity in the U band (bottom panel) and X-rays (middle panel) 
(0.3- 10 keV) with a binning time of 10 s, together with the ratio of the 
U and X-ray band (top). 


rates of this source, the number of counts in every interval is so 
low that only extremely loose 3 cr upper limits on the pulse am- 
plitude of the order of Aul ~ 0.5 could be obtained. 

In order to draw more stringent constraints on the presence 
of a coherent signal in the dataset, we applied quadratic coher- 
ence recovery techniques (QCRT; Wood et al. 1991; Vaughan 
et al. 1994). Such techniques rely on the correction of the times 
of arrival of the X-ray photons by using a quadratic time trans- 
formation, under the assumption that the sinusoidal Doppler 
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Fig. 9. Top panel: power spectrum obtained by averaging 242 intervals, 
each ^124s long, and logarithmically re-binning the spectrum with a 
factor of 1.1. The white noise level of 1.9880(4), evaluated by mod- 
elling the spectrum between 0.5 and 2 kHz with a constant, is subtracted 
from each power estimate. Bottom panel: residuals in units of cr of the 
power spectrum with respect to the best-fitting model composed of three 
lorentzians centred at zero frequency. 



Fig. 10. Right scale', upper limits at a 3<x confidence level on the pulse 
amplitude during the EPIC-pn observation of XSS J1227, as obtained 
by performing a QCRT over 1983.7 s long intervals (see text). Left 
scale', the net 0.5-10 keV light curve recorded by the EPIC-pn is also 
over-plotted. 


modulation of the arrival times introduced by the orbital mo- 
tion is well approximated by a parabola over the considered in- 
tegration time. The optimal integration time for the detection 
of a signal with QCRT was derived by Johnston & Kulkarni 
(1991) as 


( v r i/3 i 

( fj 

f m 2 \ 

- 4/15 

1 1 

( sin i \ 

V 300 Hz/ 1 

1 4.3 h/ 

(0.3 M 0 J 

(sin 45°/ 


( 2 ) 

We therefore split the data in N m u = 15 intervals composed of 
2 23 bins, each of length equal to tbi n as before, and performed 
FFT on the time series corrected with the relation t' — at 2 * * * , 
where a — -1/2 (asin//c) (27rv) (27r/P or b) 2 sin</io, asin//c is 
the projected semi-major axis of the primary star orbit and </>o 
is the orbital phase at the beginning of the time series. The 
maximum and minimum values of a, as well as the value 6a 
by which a is incremented in each of the corrections, have 
been determined by Vaughan et al. (1994) as a m ax = -<r m [ n = 
8.80 x 10" 8 (g/0.29)(P orb /4.3 hr 4/3 (M tot /1.7 M 0 ) 1/3 s" 1 , and 
6a = l/(2v Ny T 2 ). By considering the range of values for the 
parameters of the binary system and the reported uncertainties, a 
range of masses between 1.4 and 2.0 M Q for the mass of a puta- 
tive NS in this system, VNy - 21 14.35 Hz and 7 fft — 1983.7 s, 
we set a m . dx — -a m i n = 1 x 10 _7 s' 1 and 6a - 6 x 10~ u s -1 
so that N a = 3334 corrections were performed on each of the 
Mntv = 15 time intervals. The detection threshold is therefore 
P jet = 62.58 in the frequency range, where the noise is compat- 
ible with being distributed as a chi-squared with two d.o.f. (see 
Vaughan et al. 1994). No signal above this level was found. The 
upper limits at 3cr confidence level set on the pulse amplitude in 
each of the intervals are plotted in Fig. 10, together with the light 
curve in the 0.5-10 keV band. The values we find as limits on 
the pulse amplitude (~0. 15-0.25) are larger than those usually 
measured from pulsars in LMXB systems. 

3.6. Persistent X-ray and gamma-ray emissions 

The X-ray history of XSS J 1227 collected over 7 yr shows that 
it is a persistent and rather stable X-ray source, with flares and 


dips being peculiar and permanent characteristics. The associ- 
ated Fermi - LAT source 2FGL J1227 is reported in both 1st and 
2nd catalogues to be non-variable on a monthly timescale over 
the first two years of Fermi operations (see Nolan et al. 2012) 
and also on timescale >4d in the 0.1-300GeV range over the 
first ~25.5 months of Fermi-hXT observations (Hill et al. 201 1). 

We then extended the gamma-ray coverage in the 
0.1-300GeV range up to Apr. 18, 2012 and show the whole light 
curve (Fig. 11) with the same 4d binning as Hill et al. (2011). 
The average photon flux is ~1 x 1 0 7 photons cuT 2 s', slightly 
lower but consistent with that obtained by Hill et al. (2011). 
Differences are due to the different photon selection criteria and 
IRF adopted at the two epochs of data analysis 2 We have also 
overlaid the X-ray pointings since 2008, from which it is possi- 
ble to conclude that both X-ray and gamma-ray emissions are de- 
tected simultaneously. Since a shorter time bin size is precluded 
by the low photon flux, it is not possible to infer whether the 
gamma rays are emitted steadily or bursts or flares. 

4. The spectral characteristics 

In dM10 we reported the analysis of the average X-ray spec- 
trum in the 0.3-100keV range using XMM-Newton, RXTE, and 
INTEGRAL. The spectrum is featureless and well fitted with 
an absorbed power-law with T = 1.7 + 0.02. Similar results 
were found by Saitou et al. (2009). The spectral shape was 
also found to be invariant to flux variations, e.g. during quies- 
cence, quiescent dips, and flares, but not during the hard post- 
flare dips, where a partial covering absorber better described the 
spectrum (dM10). We have further checked this behaviour using 
the 201 1 XMM-Newton EPIC-pn data, fitting them with an ab- 
sorbed power-law and obtained similar results: T = 1.64 ± 0.01 
(quiescence), F = 1.65 + 0.03 (flares), F = 1.71 ± 0.04 
(quiescent dips), and F = 0.74 ± 0.08 (post-flare dips). The 
column density of the absorber is within errors the same in all 
fits and similar to that found in the 2009 XMM-Newton data, 

2 Details of reccomendation for data selection at the 

Fermi Science Support Centre: http://fermi.gsfc.nasa. 

gov/ssc/data/analysis/documentation/Cicerone/ 

Cicerone_Data_Exploration/Data_preparation.html 
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Time since MJD 54683.158 (days) 


Fig. 11. Fermi - LAT light curve in the range 0. 1-300 GeV with a 4 d bin- 
ning time from Aug. 2008 to Apr. 2012 extracted in a 1° aperture. The 
mean flux is reported together with the lcr standard deviation. The back- 
ground is not subtracted. The X-ray coverages by Suzaku (magenta), 
RXTE (green), XMM-Newton (red), and Swift (blue) are also reported. 
(See electronic version for colour figure.) 


(Vh = 0.75 ± 0.04 x 10 21 crrT 2 , and consistent with the galactic 
medium in the direction of the source. The fit to the post-flare dip 
spectrum improves when adding a partial covering absorber. We 
find ^- 2 /d.o.f. = 33/25 vs. ^- 2 /d.o.f. = 43/27 with and without 
this component. This component is found to be significant at the 
2. lcr confidence level. The column density and partial covering 
fraction of the absorber are 3.85/ 2 ^ x 10 21 cm -2 and 0.61 ± 0.2, 
respectively. Also, fits to the RXTE quiescent and flare spec- 
tra give similar power-law indexes, further corroborating the un- 
variancy of the spectrum during active and persistent level. The T" 
spectrum of the quiescent dips might indicate a softening, but _ 
it does not require an additional soft component.The 2-10keV “ 

fluxes as derived for the 20 1 1 XMM-Newton data are 0.26 + 0.04, 

3.12 + 0.1, 0.98 + 0.05, and 1.4 + 0.1 xlCL 11 ergcnD 2 s _1 during | 
the quiescent dips, the flares, the post-flare dips, and quies- 
cence, respectively. These are consistent with those found in 
2009 (dM10). 

To construct the broad-band SED we implemented the X-ray 
unabsorbed spectrum obtained in dM 10 using th e XMM-Newton, 
RXTE, and INTEGRAL data. This is reported in Fig. 12 together 
with the best-fit power-law to the quiescent X-ray average spec- 
trum. We also overlaid the extinction-corrected XMM-Newton 
OM UV and U band average fluxes (dM10), the extinction- 
corrected optical spectrum by Masetti et al. (2006), and the 
2MASS nIR fluxes (see dM10 for details on extinction). The 
UV/optical and nIR data are consistent with each other, but are 
off by more than one order of magnitude from that extrapo- 
lating the X-ray power-law. The simultaneity of the UV and 
U band data with the X-ray data then implies a different com- 
ponent contributing to the UV/optical/nlR flux. We therefore 
analysed the UV/optical/nlR SED and found that it is best fitted 
by a composite function consisting of two blackbodies at 7), = 

12 800 + 600 K and T c = 4600 ± 250 K (uncertainty at lcr c.l.) 
(y 2 /d.o.f. = 11/20). In Fig. 13 we report the UV/optical/nlR 
SED, where the optical spectrum from (Masetti et al. 2006) was 
rebinned after removal of emission lines together with the two- 
component blackbody function. The latter is also reported in 
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Fig. 12. Combined unabsorbed SED using the gamma-ray fluxes from 
Fermi-LAT (blue) and AGILE/GRID (red), the XMM-Newton (ma- 
genta), RXTE (light blue), and INTEGRAL (cyan) average X-ray spec- 
trum, the XMM-Newton OM UV and U band (black stars) fluxes from 
dM10, the optical spectrum from Masetti et al. (2006) (green), the 
2MASS nIR measures and the ATCA radio measures from Hill et al. 
(2011) (black stars). Overlaid are the simple power-law forms derived 
from the X-ray and gamma-ray spectra and from the radio-flux ratios as 
well as the best-fit composite blackbody model to the UV/optical/nlR 
fluxes. (See electronic version for colour figure.) 



Wavelength (A) 

Fig. 13. Unabsorbed UV to nIR SED, obtained combining the 
XMM-Newton OM UVM2 and U band fluxes (filled circles), the optical 
spectrum from Masetti et al. (2006) (crosses), and the 2MASS fluxes 
(empty squares) together with the best-fit composite (solid line) func- 
tion with two blackbodies with temperatures of 12 800 K and 4600 K 
(dotted lines). 


Fig. 12. The projected areas of the hot component is 5% of the 
cool one. The low-temperature blackbody could be compatible 
with that of a late-type star in the range K2-K5 V, which is of 
earlier type than what is estimated in dM10, assuming an or- 
bital period of 4.3 h. For a K2-K5V star (log g - 4.5), the 
projected area would imply a distance of 1.7 ±0.3kpc. The J 
and K band fluxes of this component are 98% and 99% the 


A89, page 11 of 16 


A&A 550, A89 (2013) 


total observed fluxes, respectively. Hence, it is reasonable to as- 
sume that the nIR flux is almost totally dominated by the donor 
star. We then also use the lowest /-band magnitude observed 
in dMIO photometry (/ = 16.9), corrected for interstellar ab- 
sorption Aj — 0.12 mag (dMIO), and the absolute magnitudes 
corresponding to a K2-K5 V star (Knigge 2006), thus obtain- 
ing a distance of 2.3-3.6kpc. If the secondary star is indeed of 
X-spectral type, it would in turn indicate a binary period in the 
range 7-9 h (Smith & Dhillon 1998). These values and the ~6 h 
variability detected in the near-UV might suggest a longer or- 
bital period. The high-temperature component instead suggests 
a non-stellar contribution from an accretion disc or flow, because 
the optical spectrum is characterized by emission lines. The pro- 
jected area of the hot emitting region is ~1 x 10 20 <7 2 kpc cm 2 . 

Assuming the association with the 2FGL source, we also 
constructed the gamma-ray portion of the spectrum using the 
AGILE/GRID upper limits and the Fermi- LAT fluxes as listed 
in the 2nd source catalogues (Fig. 12). The spectrum is also re- 
ported together with the best- fit 2FGL catalogue power-law with 
photon index 2.33. The high-energy portion of the gamma-ray 
spectrum might suggest a decay that Hill et al. (2011) fitted with 
an exponential cut-off at ~4 . 1 GeV, found to be significant at 4cr. 
The logParabola used in the 2FGL catalogue could also support 
this feature. The AGILE/GRID upper limits are roughly compat- 
ible with the Fermi-LAT measures. The gamma-ray SED does 
not support a peak in the l-100MeV range, as suggested in 
dMIO. Therefore there are no strong evidences that the X-ray 
and gamma-ray fluxes are linked. 

We furthermore include the dual band ATCA fluxes at 
5.5 GHz and 9 GHz from Hill et al. (2011), acquired in 2009, 
of the radio counterpart to XSS J 1227. Although not contempo- 
raneous, the X-ray and gamma-ray histories allow us to adopt the 
radio measures for the broad-band SED. The flux-density ratio 
at these frequencies, assuming a power-law distribution S v oc y“, 
gives a = -0.5 (Hill et al. 2011). Taken at its face value, it is 
consistent with optically thin synchrotron emission. An extrap- 
olation of this power-law to higher energies falls well below the 
X-ray spectrum but roughly matches the soft gamma-ray flux at 
~100MeV. However, we note that the large uncertainties in the 
spectral index derived by Hill et al. (2011) may also allow a flat 
radio spectrum. 


5. Discussion 

We have presented new X-ray data of XSS J 1227 and new 
gamma-ray data of its positionally associated Fermi- LAT 
counterpart to infer the long-term temporal history and the pos- 
sible nature of this intriguing source. 

5.1. A steady source: the persistent, flaring, and dipping 
states 

XSS J1227 is found at about the same flux level (/A inkcv ~ 1 x 
10~ n erg cm' 2 s~ 3 ) over 7 yr of observations, displaying pecu- 
liar short-term behaviour with flares and dips in all observations. 
Flares are simultaneously detected in both X-rays and UV band 
on timescales of -600-850 s. All X-ray flares occur in the form 
of flare-dip pairs. In dMIO we drew comparison with type-II 
bursts sources, although with different flare energetics and du- 
ration. Flares appear to occur on a timescale of <;9- 1 1 h and are 
grouped in multiple shorter events with a total duration of ^ 1 h. 
Flares were also reported to occur in the nIR band by Saitou 
et al. (2011). The large variability (~0.4mag) detected in our 


previous /-band data (dMIO) was not recognized to be in the 
form of flares due to the poor statistics. The optical photometry 
by Pretorius (2009) and dMIO also shows large-amplitude short 
term variations that could be affected by flares (dMIO). 

While X-ray flares are characterized by unvariant spec- 
tral shape with respect to the quiescent persistent level, the 
post-flare dips involve absorbing material that is dense (Ah = 
4 x 10 21 cur 2 ), partially (—60%) covering the X-ray source. The 
simultaneous occurrence of flares at lower energies in the near- 
UV corroborates the interpretation given in dMIO, where flare- 
dip pair events are due to accretion events onto the compact ob- 
ject followed by a slow refilling of the missing regions of an 
accretion disc. The emptied region should be extended because 
it is detected at UV and near-UV energies. 

On the other hand, the ordinary dips observed during quies- 
cence in X-rays and near-UV appear to occur randomly. Here we 
note that optical ground-based photometry does not reveal such 
sharp dips (Pretorius 2009, dMIO), indicating that dips origi- 
nate close to the compact star. The sharpness of their profiles, 
the different lengths (from ~ 15 s up to ~1 100 s), the timescales 
of decays/rises and the invariance (or possibly slight softening) 
of the X-ray spectrum might point to occupation due to dis- 
crete material. If this matter is located at the disc edge, the 
dip ingresses/egresses might be used to estimate the size of 
the eclipsed X-ray emitting region, similar to LMXB dippers. 
However, LMXB dippers are known to show intensified dipping 
activity at specific orbital phases with hardening of the spectrum 
due to absorbing material located close to the disc bulge (Boirin 
et al. 2005; Iaria et al. 2007; Diaz Trigo et al. 2009). The fre- 
quent and apparently non-periodic occurrence of deep dips in 
XSS J 1227 could be related to opaque material distributed ran- 
domly above the disc. Although caution should be taken, we fol- 
low Church & Balucinska-Church (2004) for LMXB dippers, 
using the longer dips where the X-ray count rate drops to zero. 
The range of ingress/egress times is ~40-90s. We assume an 
accretion disc radius Raise ~ 0.90 RLobe.i- Given the uncertainty 
in the orbital period, we use the wide range ~4— 9 h and adopt 
a mass ratio q = 0.8. With these values, the mass of the com- 
panion star would be in the range 0.4-1 M Q (Smith & Dhillon 
1998) and the disc radius would result ~4-5 x 10 10 cm. If the 
material is located at this distance and the angular size of the 
material is larger than the X-ray emitting accretion disc corona 
(ADC, see Eq. (1) of Church & Balucinska-Church 2004), the 
time of ingress is related to the radius of the emitter. In this 
case, we obtain Radc ~ 2-8 x 10 s cm. This is about one or- 
der of magnitude smaller than those determined for LMXB dip- 
pers (Church & Balucinska-Church 2004). The size of the ADC 
scales with the X-ray luminosity (Church & Balucinska-Church 
2004) would imply Lx ~ 0. 5-1.0 x lO^ergs^ 1 . This is about 
two orders of magnitude larger than that derived for a distance 
of 1 kpc (dMIO). Even at distances of 2.4-3.6kpc, the X-ray 
(bolometric) luminosity is lower (3 - 8 x 10 34 erg s _1 ). Therefore, 
either the disc corona is limited in extension or the occulting 
material is not located at the adopted disc rim. 

In addition, the dips are observed in the UV, implying that 
the hot 13kK component that dominates the UV light also 
suffers occultation. This component has an emitting area of 
7t(0.2R o ) 2 Z 2 kpc - This is smaller than the adopted disc radius 
but larger than those determined for similar hot components in 
LMXBs and attributed to the bright disc hot spots or a region 
closer to the centre (Froning et al. 2011; Hynes & Robinson 
2012). The UV/optical emitting region is larger than that esti- 
mated for the X-ray emitter. This can explain why the UV flux 
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deficiency is lower (~40%) than that in the X-ray band during 
the dips. 

Furthermore, the fact that the ratio of bolometric fluxes be- 
tween the X-ray and UV components is ~6 and the X-ray peak 
intensity at flares is about twice that in the UV indicates strongly 
that the UV light originates from X-ray reprocessing in a larger 
region surrounding the X-ray source. 

In summary, the long-term X-ray history of XSS J 1227 over 
7 yr and the gamma-ray light curve of the associated source 
2FGL J1227 over 4.7 yr indicate that this peculiar binary is a per- 
sistent source in both energy ranges and that it always displays 
the same X-ray variability characteristics. We cannot, however, 
exclude the possibility that the gamma rays are also emitted in 
the form of flares with shorter timescales than 4 d. 

5.2. An LMXB hosting an MSP? 

Analogies to the MSP binary PSRJ 1023 +0038 (Tam et al. 2010) 
were suggested by Hill et al. (2011) because of similarities of 
the optical spectra (both sources were formerly proposed as 
CVs), the existence of radio and Fermi-LXT counterparts and 
the short, putative ~4 h orbital period, and large optical flicker- 
ing. PSR J1023+0038 was proposed to be active in gamma rays 
when accretion switches off and evidence of accretion stopping 
at least once comes from an optical, although non-simultaneous, 
spectrum (see Tam et al. 2010 and reference therein). The pos- 
sibility that accretion switches off in XSS J 1227 was left open 
in Hill et al. (201 1), but the present study strongly indicates that 
the source has always been in an accretion state. This is further 
confirmed by recently acquired optical spectroscopy (de Martino 
et al., in prep.). 

We draw a further comparison with MSP binaries detected 
by Fermi. The gamma-ray spectra of PSRJ1023+0038 and 
2FGLJ1227 are rather similar, with a steep power-law in- 
dex F2.5 (in PSR J 1023 +003 8) and F = 2.3 in 2FGLJ1227. 
Although not strongly constrained, a power-law exponential cut- 
off model for PSR J 1023+0038 gave F = 1.9 and £ cutoff = 
0.7 GeV (Tam et al. 2010) as compared with F = 2.2 and 
Cutoff = 4.1 GeV inferred for 2FGL J1227 by Hill et al. (201 1). 
Ransom et al. (2011) fitted the gamma-ray spectra of PSR J0614- 
3329, PSRJ1231-1411 and PSR J22 14+3000, with a power-law 
exponential cutoff, obtaining £ C utoff in the range 2. 5-4. 5 GeV 
and F ~ 1.4. Kerr et al. (2012) find for PSR J0101 -6422 E catoS = 
1.9 GeV and T = 0.9. Also for the very recently discovered 
black- widow MSP PSR J18 16+45 10, Kaplan et al. (2012) found 
T = 2 and £cutoff ~ 7.5 GeV. Therefore, based on the gamma-ray 
spectrum, 2FGL J1227 is not much different from the variety of 
MSP binaries recently discovered by Fermi. 

The X-ray (0.2-100 keV) to gamma-ray (0.1-100GeV) lu- 
minosity ratio of our source is ~0.8, whilst most of the MSP 
in LMXBs have lower ratios by 2-3 orders of magnitudes 
(Bogdanov et al. 2006, 2011; Ransom et al. 2011). Our explo- 
ration of the overall SED also indicates that the X-ray flux, 
dominated by a non-thermal component, is higher than that ex- 
pected from rotation-power MSP emission. In the latter, the ther- 
mal flux from the NS atmosphere either dominates (Ransom 
et al. 2011) or is detectable, as in the case of PSR J 1023+0038 
with a thermal fraction ~3-6% (Bogdanov et al. 2011; see 
also Takata et al. 2012). The non-thermal X-ray emission, if 
originating from an intrabinary shock, such as in the black- 
widow MSPs, would be the extension of the gamma-ray flux at 
these energies (Takata et al. 2012), and this is not the case for 
XSS J1227. Intrabinary shock as source of non-thermal emission 


in the quiescent MSP binary PSR J1023+0038 was proposed by 
Bogdanov et al. (201 1). 

Our search for X-ray msec pulses in XSS J1227 only pro- 
vided upper limits ~15% to pulse amplitudes. In accreting 
MSPs, pulse amplitudes in excess of 10% are rarely ob- 
served, but fractional amplitudes can range between 0.43% and 
~30-40% (see review by Patruno & Watts 2012). Here we note 
that in PSR J1023+0038 the X-ray pulse amplitude was found 
to be 1 1 + 2% (Archibald et al. 2010). Therefore, we cannot ex- 
clude that XSS J 1227 harbours a fast-spinning pulsar. However, 
in non-accreting MPS the X-ray pulses would arise from the NS 
atmosphere or hot spots that we have not detected in the X-ray 
spectrum of XSS J 1227. A non-thermal component arising from 
intrabinary shock would only produce an X-ray orbital modu- 
lation, as is the case of PSR J1023+0038, which shows large 
amplitude orbital variability (Bogdanov et al. 2011). While in 
the 2009 XMM-Newton data, a possible weak X-ray variability 
could be present, this is not detected in the higher quality EPIC- 
pn data in 2011. Hence, all this points to an accretion-driven 
X-ray emission. 

Accretion acting in this source is further corroborated by the 
UV/optical/nIR portion of the SED, which reveals two spectral 
components: a hot one at ~ 1 3 kK and a cool one at 4.6 kK. A hot, 
but single, spectral component was also recently identified in the 
UV/optical SED of the MSP binary PSR J1 8 16+45 10 (Kaplan 
et al. 2012), making it the black-widow binary pulsar hosting 
the hottest low-mass companion (^ 1 0 kK). We however note 
that no nIR data are available to date for PSR J 1 816+4510. In 
XSS J1227, the cool component is instead revealed from the nIR 
2MASS measures and its temperature suggests the contribution 
from a late-type K2-K5 companion, earlier than that estimated 
in dM10. Hence, the UV/optical component can be reasonably 
ascribed to accretion. The temperature is, however, higher than 
that of accretion discs (see Hynes & Robinson 2012) and the 
size of the emitting region could suggest that it originates closer 
to the compact object. The X-ray and UV bolometric luminosity 
ratio is ~6, compatible with an X-ray heated region. We note that 
this ratio has been proposed to be a discriminant between NS and 
BH binaries, with NS having X-ray-to-UV flux ratios about ten 
times higher than BH systems (Hynes & Robinson 2012). 

If the radio spectral power-law index is truly a = -0.5, then 
the extrapolation to the gamma-ray regime roughly matches the 
observed flux at ~100MeV. However, searches for fast pulses 
in the radio counterpart of XSSJ1227 failed (Hill et al. 2011). 
Therefore, while the flux ratio could be consistent with non- 
thermal emission, the lack of detection of pulses does not favour 
a rotation-powered radio pulsar. 

5.3. An LMXB powering a jet? 

The hard X-ray spectrum and the low-luminosity of XSS J1227 
indicate an LMXB in a deep faint hard state. From the detec- 
tion of quasi-simultaneous X-ray and nIR flares and the low 
X-ray luminosity, Saitou et al. (2011) argued that XSSJ1227 
could be a microquasar with a synchrotron jet in a prolonged 
low-luminosity state at 1 0 4 Lem for a stellar mass BH or NS at 
1 kpc. As a consequence of synchrotron cooling in an optically 
thin medium, the energy break expected in these binaries is at 
~500keV (Saitou et al. 201 1), much lower than that detected at 
~4GeV in 2FGLJ1227. 

High-energy gamma rays were detected in galactic bi- 
naries hosting high-mass early-type stars, such as LS 5039, 
LSI +61°303 (Abdo et al. 2009a, b; Hadasch et al. 2012), 
1FGLJ101 8.6+5856 (Ackermann et al. 2012), and CygX-3 
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during flares (Tavani et al. 2009b; Abdo et al. 2009c). The latter 
is the only gamma-ray source with a BH candidate while the oth- 
ers are proposed to possibly harbour NSs. The gamma-ray spec- 
trum of LSI+61°303 is compatible with a power-law exponen- 
tial cut-off at ~4GeV (Abdo et al. 2009a; Hadasch et al. 2012), 
while LS 5039 shows a power-law exponential cut-off at ~2 GeV 
(Abdo et al. 2009b; Hadasch et al. 2012). The nature of GeV 
emission from these binaries is still unclear and controversial 
(Hadasch et al. 2012; Papitto et al. 2012; Torres et al. 2012). 
They could be powered either by magnetospheric emission, as 
seen in many gamma-ray detected pulsars, or by electrons in the 
jets accelerated up to TeV energies (leptonic model; Hadasch 
et al. 2012). For CygX-3, the Fermi flux is compatible with the 
extrapolation of the X-ray tail to 100 MeV and hence is believed 
to be the high-energy tail of the soft gamma-ray spectral compo- 
nent. Also, in this case the origin of gamma-ray photons is un- 
clear and is proposed to be leptonic (Abdo et al. 2009c). Since no 
LMXB with a radio jet has been detected in the gamma rays so 
far, it is difficult to predict the origin of high-energy gamma-ray 
emission. 

We alternatively explored the radio-to-X-ray emission from 
LMXBs hosting NS and BH in the hard faint quiescent states. 
The persistent radio emission from BH binaries during these 
states is generally consistent with a flat or slightly inverted 
spectrum that, in analogy with compact extragalactic radio 
sources, is attributed to the superimposition of a number of 
peaked synchrotron spectra generated in a jet (Gallo 2010). Atoll 
NS binaries in low-hard states (island states) have also been de- 
tected at radio frequencies, though most are upper limits with 
a handful of systems detected simultaneously with the X-rays 
(Migliari & Fender 2006). Usually, dual-band radio spectra are 
not well constrained and could be consistent with flat (a - 0) or 
optically thin a ~ -0.6 slopes, such as Aql X-l or 4U 1728-34 
(Migliari et al. 2003; Migliari & Fender 2006). In NS binaries, 
jets can nevertheless be powered during these low states but are 
about one order of magnitude less powerful than BH jets. The ra- 
dio spectrum ofXSSJ 1227 is not well constrained, a = 0.5 + 0.6 
(Hill et al. 2011), leaving the possibility of a similar origin as 
those in BH or NS binaries. 

The BH binaries spend most of their time in such states, 
which could be compared to the prolonged low-state of 
XSS J1227. The X-ray emission from these binaries originates in 
an optically thick accretion disc with a hard component, which 
could be an ADAF-like (advection dominated accretion flow; 
Narayan & Yi 1994) or from a jet. The jets are coupled with 
the accretion flow, although the disc-jet connection is still un- 
clear. The radio emission from the jet was found to be corre- 
lated with that in the X-rays during the hard state: F ra dio °c /•’*, 
with b ~ 0.5-0. 7 (Gallo et al. 2003). A similar relation was 
found for atoll NS binaries in the island states but with a steeper 
correlation coefficient b ~ 1.4 (Migliari & Fender 2006). Hence, 
the X-ray emission in the BH systems is radiatively inefficient, 
with an approximate relation Lx °c M 2 , whilst the jet power in 
both classes scales linearly with the accretion rate. Keeping in 
mind the large uncertainty in the radio spectral slope, we have 
explored the possibility that XSS J 1227 could host a compact 
jet following Coriat et al. (2011), where the radio/X-ray lumi- 
nosities of BH binaries and NS LMXBs are compared. This is 
depicted in Fig. 14, adapted from Fig. 5 of Coriat et al. (2011). 
Here, the radio (9 GHz) luminosity against the X-ray 3-9 keV 
unabsorbed luminosity is shown for a sample of BH candidates 
in LMXBs (GX 339-4, H 1743-322 and V404Cyg) during hard 
X-ray states and for the two atoll NS binaries in the island state. 



Fig. 14. Radio luminosity at ~9GHz versus 3-9 keV X-ray luminos- 
ity for the BH candidates GX 339-4, H 1743-322 and V404Cyg dur- 
ing their hard X-ray states, together with the atoll-type NS binaries 
Aql X-l and 4U 1728-34 during their island states (see footnote 2 for 
details on the data). The radio/X-ray luminosity correlation F ra di 0 °c 
with b ~ 0.6 for BH binaries and b ~ 1.38 for island state NS bina- 
ries is also reported. The positions of XSS J 1227 at the three distances 
1 kpc, 2.4 kpc, and 3.6 kpc are marked in red. Adapted from Coriat et al. 
(2011). (See electronic version for colour figure.) 


Aql X-l and 4U 1728-34 3 In the radio vs. X-ray luminosity 
diagram, the BH candidate H 1743-322 moves to lower X-ray 
luminosities, transiting the two correlations and joining the stan- 
dard BH correlation at low luminosities (Coriat et al. 2011). This 
is the first BH candidate to follow both correlations depending 
on the X-ray luminosity. 

We locate XSSJ1227 in this diagram by using the ATCA 
9 GHz and 3-9 keV fluxes adopting the three values for the dis- 
tance, 1 kpc, 2.4 kpc, and 3.6 kpc. For the sake of caution, these 
are marked as upper limits because of the uncertainties in the 
radio slope and hence in the contribution of a jet. XSS J 1227 ap- 
pears to be rather away from the NS binaries locus and, in case of 
a jet in this system, it is closer to the BH correlation. Its position 
could move upwards, matching at 10 kpc the transition region 
of H 1743-322, where the latter joins the standard BH correla- 
tion. While we cannot infer the true nature of the compact star 
in XSS J1227 based on this diagram alone, if a jet is present, 
it could possibly favour a BH LMXB. If so and similarly for 
H 1743-322 (Coriat et al. 2011), this binary would be in an ex- 
tremely deep low hard state, where the accretion flow becomes 
radiatively inefficient. This flow is then dominated by advection 
in which a significant fraction of the energy is advected and not 

3 The radio-vs.-X-ray luminosity diagram is constructed using data for 
GX 339-4 from Corbel et al. (2003) and a distance of 8 kpc, for H 1743- 
322 from Table 1 of Coriat et al. (2011) and distance of 8 kpc, and for 
V404 Cyg from Corbel et al. (2008) and a distance 2.39 kpc. For the 
atoll NS binary Aql X-l, we use data from Miller-Jones et al. (2010) 
and distance 5 kpc, and for 4U 1728-34 we use data from Migliari & 
Fender (2006) and distance 4.6 kpc. 
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radiated away (Narayan & Yi 1994). The energy either crosses 
the BH event horizon or is expelled in outflows. Systems dom- 
inated by jet emission, where most of the energy is channeled 
in a jet, emit X-rays at the base of the jet. These are produced 
either by synchrotron emission or by inverse Compton scatter- 
ing by outflowing particles. In the previous section, we derived 
an approximate estimate of the size of the X-ray emitting region, 
which is much smaller than that of typical ADC in LMXBs. 
Coriat et al. (2011) also suggest that H 1743-322 possesses two 
components, one radiatively efficient and one inefficient, which 
co-exist and dominate alternatively in the X-ray band during the 
hard state. In the case of XSS J1227, the bulk of X-rays does not 
seem to originate in a jet because of the evidence of accretion. 
Further radio data would be important to assess the true nature 
of the radio emission. 

6. Conclusions 

We have presented new X-ray data to study the long-term history 
of XSS J1227 and to infer whether the X-ray emission is variable 
on timescales down to msec in order to shed light on the nature of 
this peculiar source. We here summarize the main conclusions: 

- The long-term X-ray history of XSS J 1227 over a time- 
span of 7 yr shows that it is a persistent low X-ray lu- 
minosity source and not a transient one. The associated 
source 2FGL J1227 is also found to be a steady high-energy 
gamma-ray source over an overlapping period of 4.7 yr. 
The co-existence of X-ray and gamma-ray emissions does 
not favour an MSP binary nature if the two sources are 
associated. 

- XSSJ1227 displays flares and dips in all of the X-ray ob- 
servations, with flares possibly occurring on a timescale of 
9- 1 1 h and duration of ~ 1 h. Flares are grouped in multi- 
ple events. We confirm previous finding of the occurrence of 
flare-dip pairs. The spectral hardening in these dips are due 
to an intervening dense, Nn ~4x 10 21 cur 2 , absorber cov- 
ering ~60% of the X-ray source. Intense dipping activity is 
found during quiescence. These dips have variable lengths 
from ~10s to 1100 s. A marginal spectral softening is de- 
tected in the longer dips. We estimate the size of the X-ray 
emitting region to be ~2-8x 1 0 X cm, much smaller than those 
derived for LMXB ADC dippers. The UV flares and dips are 
found to occur simultaneously with the X-ray ones. 

- We searched for fast pulsations down to msec in the fast- 
timing X-ray data obtaining upper limits to fractional pulse 
amplitudes of 15-25%. These do not exclude an MSP in this 
system. A rotation-powered MSP could be favoured by the 
SED at gamma-ray energies and possibly at radio frequen- 
cies, but not from the X-ray and UV/optical/nIR data. The 
X-ray source is consistent with being an accretion-powered 
binary. 

- From the combined UV/optical/nIR spectrum, we find ev- 
idence of a hot component at ~13kK and a cool one at 
' 4.6 kK. The former is consistent with an X-ray-heated ac- 
cretion region, which is smaller than that of a tidally trun- 
cated accretion disc but larger than those found in other 
LMXBs with either NS or BH. The cool component would 
suggest a late-type K2-K5 companion star and in turn a dis- 
tance of 2.4-3. 6 kpc. Then the orbital period would be in 
the range of 7-9 h. This would embrace the 4.3 h periodicity 
previously inferred (dMIO) if it were the first harmonic. A 
20% variability at 6.4 h is found in the quiescent near-UV 
photometry, further suggesting a long orbital period. 


- The X-ray-to-UV luminosity ratio is compatible with an 
LMXB hosting an NS that is much larger than that observed 
in BH LMXBs. On the other hand, the uncertainties in the 
radio spectral slope may also allow the possibility of a com- 
pact jet. In this case, the radio and X-ray luminosities might 
indicate XSS J1 227 as a BH binary rather than an NS atoll 
system during island states. The persistent low luminosity of 
this source could be due to a radiatively inefficient accretion 
flow, where the radio emission originates in a jet but the 
X-ray emission does not. If a BH binary, XSS J1227 were the 
first LMXB to be associated with a high-energy gamma-ray 
source. 

To further progress on the intriguing nature of this source, it is 
crucial both to obtain new observations in order to infer whether 
flares also occur in the radio domain and to study the radio spec- 
trum and optical spectroscopy in order to determine the true 
orbital period. 
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